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This paper describes the development of an experiment to acquire data for developing 
and validating computational fluid dynamics (CFD) models for turbulence in supersonic 
combusting flows. The intent is that the flow field would be simple yet relevant to flows 
within hypersonic air-breathing engine combustors undergoing testing in vitiated-air 
ground-testing facilities. Specifically, it describes development of laboratory-scale hardware 
to produce a supersonic combusting coaxial jet, discusses design calculations, operability and 
types of flames observed. These flames are studied using the dual-pump coherent anti- 
Stokes Raman spectroscopy (CARS) - interferometric Rayleigh scattering (IRS) technique. 
This technique simultaneously and instantaneously measures temperature, composition, and 
velocity in the flow, from which many of the important turbulence statistics can be found. 
Some preliminary CARS data are presented. 


I. Nomenclature 

Variables: 

ER = Equivalence ratio 

H = Convective heat transfer coefficient 

M = Mach number at nozzle exit 

M c = Convective Mach number 

Menthaipy = Flight Mach number with the same enthalpy 

q w = Surface heat flux 

r = Radial distance 

T= Nozzle exit static temperature 

T wa u = Burner interior wall temperature 

x = Axial distance 


Flame state labels: 

NF = No flame 

DF = Detached flame 

EAF = External attached flame 

BAF = Base attached flame 

M = Flame marginal stability 

0 = Burner would not light 


II. Introduction 

C omputational fluid dynamics (CFD) methods that employ the Reynolds-averaged Navier-Stokes (RANS) 
equations are widely used in the design and analysis of hypersonic airbreathing engine flow paths. These 
methods require models of statistical quantities of the turbulence, which have to be empirically calibrated and 
validated. In particular, new models for turbulent Schmidt and Prandtl number, as well as for turbulence-chemistry 
interactions, are required. 1 While suitable data is available for low-speed flows with combustion, data is still lacking 
in supersonic combustion. Goyne et al. 2 report measurements using particle-imaging velocimetry of mean 
streamwise velocity in a dual-mode scramjet. At NASA Fangley, several data sets have been acquired in a H 2 
fueled supersonic combustor using the coherent anti- Stokes Raman spectroscopy (CARS) technique 3 and the dual- 
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pump CARS technique originally developed by Robert Lucht and coworkers 4 ’ 5 and extended by us. 6,7 The former 
technique gave temperature only while the latter gave both temperature and composition. Data included both mean 
flow and turbulent statistics, although the uncertainty in the latter was high due to both instrument precision and the 
number of measurements from which the statistics were formed. International work in this area includes 
measurements in scramjet combustors conducted at ONERA (France) and DLR (Germany) using CARS 8 , and other 
non-intrusive techniques. 

Available data sets are limited to only a subset of the important variables (which are temperature, composition, 
and velocity) and data sets that include accurate turbulence statistics do not exist. To meet the need for 
instantaneous velocity measurements to complement dual-pump CARS measurements of temperature and 
composition, the interferometric Rayleigh scattering technique is being developed. 9,10 This technique uses Rayleigh 
scattering from one of the CARS laser beams to measure velocity in the same instant as the CARS measurement. 
Analysis of streams of such instantaneous measurements allows us to form the statistical quantities (means, 
variances, covariances) required by the modelers. 

Experimental facilities to provide suitable flow fields are also being developed. 11,12 An axisymmetric coaxial 
free jet has been selected since it can provide the good optical 
access required for the Rayleigh technique, and symmetry can 
be taken advantage of to minimize the number of spatial points 
required. (In order to form accurate statistics we need large 
numbers of measurements at each given point.) The flow field 
needs to provide data that is relevant to both H 2 - and 
hydrocarbon-fueled hypersonic air-breathing engines and the 
testing of such engines in ground facilities that employ vitiated 
air (the products of combustion of either H 2 or a hydrocarbon, 
enriched with 0 2 to the same content as in air). The facility 
needs to provide supersonic flows to establish the effects of 
compressibility. Finally, it is desirable to be able to create both 
supersonic combusting flows in which the flame is attached to 
the burner (flame held) and flows in which the flame is 
detached, since both types of flow are observed in hypersonic 
engine combustors. 

Two different sizes of experimental hardware are being 
developed. The first is a laboratory-scale supersonic burner, 
with flow rates compatible with testing in a combustion 
diagnostics laboratory. The center jet nozzle exit diameter of 
this device is 10 mm. This device is being used to develop the 
test matrix of experiments and to verify the CARS -IRS 
techniques. However the CARS measurement volume is -1 .5 
mm long, of the same order of magnitude of some of the length 
scales of the flow. Thus, the capstone experiments of this 
project will be performed in NASA Langley’s Direct Connect 
Supersonic Combustion Test Facility (DCSCTF) with hardware 
that is scaled up by a factor of 6.35 from the smaller device in 
its important flow dimensions (although the mechanical design 
is different). The present paper describes the development of 
the experimental hardware for the laboratory device, the test 
matrix used and the verification of the CARS-IRS 
instrumentation. 
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Figure 1. Sectional view of laboratory 
supersonic burner. 


III. Supersonic Burner 


A. Description 

Figure 1 shows the laboratory-scale burner, sectioned along the axis, shown without bolts, gas supply lines, 
spark plug and other fittings; Figure 2 is an image obtained obtained near the nozzle exit during operation, showing 
a Mach 2 jet of vitiated air and the laser beams of the CARS-Rayleigh system. The burner consists of a combustion 
chamber, 25.4 mm internal diameter and 152 mm in length, formed by an exterior “jacket” and an interior “shell” 
between which cooling water flows. An axisymmetric piece contains a central nozzle with an exit diameter of 10 
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mm. This “nozzle” is interchangeable and convergent (M <1), or supersonic convergent-divergent (M= 1.6 and 
M= 2) nozzles are available. Nozzles are designed using the axisymmetric method of characteristics with fixed ratio 
of specific heat capacities. 13 The nozzle is trapped by an “annulus” and an annular coflow nozzle passage is formed 
between it and the annulus. The coflow nozzle passage is convergent with exit width (in the radial direction) of 1 
mm. The annular base region of the nozzle, between the central jet and coflow, is normal to the axis and 3 -mm 
wide. Reactants at ambient temperature are delivered to the burner by the “injector”, which contains a central tube 
through which gaseous fuel flows, while a mixture of 0 2 (or sometimes N 2 ) and air flows in a concentric passage. 
This coflow of reactants enters the burner through a sudden expansion. The injector, jacket, shell, and annulus are 
fabricated from AISI 310 stainless steel and the nozzle is 
fabricated from Hexaloy grade SA silicon carbide. 14 

Various combinations of reactants are possible. In one 
set of experiments H 2 or C 2 H 4 fuel, air, and 0 2 are reacted 
to provide vitiated air at various temperatures (dependent 
on flow rates). The coflow is of unheated H 2 or C 2 H 4 . In 
another set of experiments the reactants to the burner are 
H 2 , 0 2 (sometimes N 2 is used instead of 0 2 ), and air in 
such a ratio that the products contain excess H 2 (rather than 
excess 0 2 ). The coflow in this case is air. In both cases, 
the resulting coaxial jet flow will mix and may react 
depending on nozzle exit Mach number and temperature. 

If reaction does take place, then the flame may be held at 
the burner or stand off from it, depending again on 
temperature and Mach number. Cases for the study of 
mixing without reaction in the jet are also possible by 
employing vitiated air, or products containing reduced 
excess 0 2 , both with air coflow. 

B. Test Matrix Calculations 

Burner stagnation and nozzle exit conditions are calculated from flow rates, heater pressure, and nozzle 
minimum and exit areas using an inviscid quasi-one-dimensional (ID) analysis 15 . The flow in the burner is assumed 
to be fully mixed (uniform in composition and thermodynamic state) and in thermodynamic equilibrium before 
entering the nozzle. The heat loss to the burner structure and cooling water is assumed, based on experience in the 
DCSCTF, to be 20% of the sensible enthalpy increase (i.e., the heat that would be released in bringing the products 
at constant pressure back to the temperature of the entering reactants). Thus, the product’s enthalpy is reduced by 
this amount. The nozzle flow is computed by ID analysis, assuming isentropic, choked flow of frozen composition 
in the nozzle. Nozzle exit conditions are computed from the geometrical nozzle throat and exit area and neglect any 
displacement effect of boundary layers. Nozzle exit pressure is forced to be one atmosphere and reactant 
temperature is assumed to be 298 K. All significant minor species are included. Other calculations assuming 
equilibrium rather than frozen chemistry in the nozzle result in significant differences in minor species composition, 
but insignificant differences (<1%) in major species, temperature, pressure, and velocity. The nozzle exit 
thermodynamic conditions for vitiated air products are related to an equivalent flight Mach number. The equivalent 
flight Mach number has the same total sensible enthalpy (sensible enthalpy based on static conditions relative to a 
standard state plus kinetic energy), where flight is assumed to be in the atmosphere on a constant dynamic pressure 
trajectory of 1000 pounds force per square foot. 

“Convective” Mach number (M c ) has been computed at the nozzle exit by the method of Papamoschou and 
Roshko. 16 These authors were able to correlate the growth rate of planar two-stream mixing layers, when normalized 
by the growth rate of incompressible layers of the same velocity and density ratio, in terms of M c . Two convective 
Mach numbers were defined, being the magnitude of the difference between the convective velocity of the mixing 
layer coherent structures and one of the streams, divided by the speed of sound of the respective free stream. A 
simple analytical method gave convective Mach numbers for each stream that were nearly equal. When correlated 
with M c , growth rates normalized by the growth rate of the incompressible layers were found to fall from 1 .0 for M c 
< 0.2 to ~0.2 for M c > 0.7. In our computations of M c , the center jet is assumed to be the high speed stream and the 
(stationary) ambient air the low. Our axisymmetric, annular flow does not closely resemble the simple planar two- 
stream mixing layer; M c here is used simply as a guide to the importance of compressibility. 
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Figure 3 Computed conditions at exit of center jet nozzle for various cases of vitiated air, with M of 0.75, 1, 
1.6 and 2, and with vitiate produced by combustion of H 2 and C 2 H 4 ; plotted as a function of equivalent flight 
Mach number, (a) Static temperature, (b) velocity, (c) M c ; mole fractions respectively of (d) N 2 , (e) H 2 0 and 

(0 co 2 . 

A matrix of calculations has been performed for cases in which the center jet is vitiated air (that is contains 23% 
by mass 0 2 ). The burner fuel is either H 2 , CH 4 (methane), or C 2 H 4 (ethylene). Methane was later dismissed as a 
practical fuel since it was found not to react in our burner, at least over the range of reactant compositions we were 
interested in studying. Ethylene has the disadvantage that, although gaseous and more reactive than methane, has a 
relatively low vapor pressure and under certain circumstances can form droplets in the supply lines. Figure 3 shows 
the results for nozzle-exit conditions. Figure 3(a) shows the nozzle-exit temperature (7), which is seen to increase 
rapidly with increasing equivalent flight Mach number, and decrease with nozzle Mach number. There is also some 
effect of burner fuel, with temperatures increasing with carbon content. Figure 3(b) shows the nozzle-exit velocity, 
which is seen to increase rapidly with both equivalent flight Mach number and nozzle Mach number. There is little 
effect of burner fuel on velocity. Figure 3(c) shows that M c depends principally on nozzle Mach number, being in 
the fully compressible range at M=\ .6 and 2, and at the upper end of the transitional range for M=l. Molecular 
nitrogen mole fraction (Fig. 3d) is dependent primarily on equivalent flight Mach number, while H 2 0 (Fig. 3e) and 
C0 2 (Fig. 3f) mole fractions depend on equivalent flight Mach number and carbon content of the fuel. 
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A similar matrix of calculations has been performed for cases in which the center jet contains excess H 2 (rather 
than 0 2 ) in mole fractions of 33%, 50% and 67%, the so called “vitiated hydrogen” cases. The burner fuel is, of 
course, H 2 . Figure 4 shows the results for nozzle exit conditions. Nozzle exit velocity (Fig. 4a) increases with T, M, 
and H 2 content, and is generally higher than with vitiated air. Convective Mach number is again primarily 
dependent upon M (Fig. 4b), being in the fully compressible range for the M= 1.6 and 2. Water vapor content 
increases with T and M (Fig. 4c). 
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Figure 4 Computed conditions at exit of center jet nozzle for various cases of vitiated hydrogen, with M of 
0.75, 1, 1.6 and 2; plotted as a function of static temperature at nozzle exit, (a) Velocity, (b) M c , (c) mole 
fraction of H 2 0. 


As has been seen, a wide range of variation of parameters of the center jet is possible. Additionally, the coflow 
composition and rate of flow may be varied. An overall equivalence ratio, ER, is defined for cases with vitiated air 
as the rate of coflow fuel to exactly consume the excess 0 2 in the center jet. Typically, ER is 0.5 or 1.0. 
Experimentally, it was found (based on IR imaging) that cases with CH 4 coflow would not react significantly for our 
test conditions, although visible light emission suggested that some dissociation of the fuel was taking place. If the 
coflow was air, the rate of flow of air was computed so as to provide sonic flow at the coflow nozzle exit with 
pressure matched to the atmosphere. This latter condition was believed, based on CFD analysis, 11 to be conducive 
to flame holding at the base for those experiments in which the center jet was vitiated hydrogen. 


C. Detailed Design Calculations 

Detailed design calculations of the burner were performed for a high heat transfer rate case: T = 1200 K, M= 2, and 
50% unreacted H 2 . The objectives of these calculations were to determine whether combustion would stabilize 
within the burner, whether the flow at the exit of the burner was uniform, and whether the burner component 
temperatures and stresses exceeded material limits. Two types of calculation were performed: calculations of the 
internal flow of the burner using computational fluid dynamics, and calculations of the temperature and stresses in 
the burner structure using finite elements analysis. These calculations were coupled via a computed surface heat 
transfer coefficient distribution. 

The Vulcan CFD code 17 was used to solve the Favre-averaged Navier-Stokes equations on a structured grid 
using the finite-volume method. The flow is assumed to be axisymmetric, and mixtures of thermally perfect gases 
were assumed. Inviscid fluxes were calculated using the Kappa = l/3 rd MUSCL scheme with the Low Diffusion 
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Flux Split Scheme of Edwards 18 , while viscous fluxes were evaluated using 2 nd order central differences. A 
diagonalized approximate factorization scheme was used for iterating the unsteady equations in pseudo-time to a 
steady-state solution. The grid, consisting of 3 1 ,248 grid points, was generated by a separate, commercial code. 

The gridded volume terminated at the nozzle exit. Grid points were clustered near the walls of the nozzles to resolve 
the boundary layers, as appropriate for the use of wall functions, and in the throat and expanding region of the 
nozzle. Turbulence was modeled using the Wilcox k-co model (1998). 19 Calculations were performed at two 
isothermal surface conditions: the first with a constant 900 K throughout, the second with 900 K for the steel 
components switching to 1650K for the silicon carbide nozzle. Reactants were assumed to enter the burner in a 
premixed state as sonic flow through a 0.1 inch diameter hole on the axis. Entering reactant (static) properties were 
1.463xl0 6 Pa, 248 K, and mass fractions of H 2 , 0 2 , and air respectively of 0.130, 0.227, and 0.643, computed as 
described in Section IIIB. Total pressure and temperature conditions were specified at subsonic inflow/outflow 
planes, while the code switched to extrapolation where the outflow is supersonic. 




Figure 5 CFD calculated properties at the nozzle exit for the nominal J=1200 K, M=2 case with T wat f=9 00 K; 
properties are plotted on the horizontal axis as a function of radial distance on the vertical, (a) M, (b) T. 

As the burner is currently (experimentally) configured, the reactants are not premixed when they enter the burner 
(see Fig. 1 and Section IIIA), as assumed in the computations. The computed configuration was in fact attempted 
experimentally. However, flashback into the supply lines occurred during the startup transient due to the pressure in 
the burner greatly exceeding the steady state immediately after igniting the reactants in the burner. It was possible 
to work around this problem by, for example, lighting at low fuel-air equivalence ratio and then increasing flows to 
the final condition. And a premixed flow entering the burner was desirable because it simplified computations and 
guaranteed uniform composition of the products. However, the desire for robust and safe operation outweighed 
these other considerations. 

Some of the results of our calculations are shown in Figs. 5 and 6. Figure 5 shows properties at the nozzle exit. 
The Mach number (Fig. 5a) is close to the design value of 2.0 in the freestream, and falls to zero at the wall, as 
expected. The static temperature (Fig. 5b) is a little higher than 1200 K at the axis and drops to 900 K at the wall. 
The 1200 K design value represents the mass flow weighted mean (accounting for heat losses to the walls of the 
burner), so the high value at the axis for this computation is to be expected. The calculated nozzle exit pressure 
varies between 107.12 kPa and 107.34 kPa, a small range but higher than the design 101.3 kPa. This discrepancy is 
attributed to the effect of the boundary layer displacement, which reduces the effective flow area at the nozzle 
throat. 

Figure 6 shows the surface heat transfer coefficient, plotted as a function of axial distance from the nozzle exit, 
computed based on the surface heat flux, q w , and an assumed recovery temperature of 2000 K: . 

H = qJ(2000-T w ) (1) 

The results show a tolerable collapse of the two different wall temperatures cases. As our later calculations will 
show, the actual nozzle wall temperature lies between 900 K and 1650 K, and 900 K is a reasonable average value 
for the burner. Thus, the computed H(x ) for the T wa if= 900 K case forms a reasonable conservative estimate for the 
calculations that follow. 
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The temperature and thermal stress 
distributions in the burner assembly were 
computed using a commercial CAD solid 
modeling and finite-elements package that also 
included a flow solver. A number of calculations 
were performed as the hardware design was 
iterated. One significant issue was the manner in 
which the silicon carbide nozzle was held. If the 
contact surface between the steel components of 
the burner and the nozzle was rigid (“bonded”) 
and conducting, stresses were very high in the 
nozzle (greatly exceeding the limits of the 
material). This had two causes. First, severe 
temperature gradients in the shell caused the 
contact surface with the nozzle to tend to distort. 

Second, conduction from the nozzle induced 
temperature gradients within the nozzle itself, 
inducing thermal stresses. These problems, once 
identified, were easily mitigated by using soft, 
graphite gaskets which, in addition to providing a 
seal, absorb the distortion and somewhat increase 
thermal contact resistance. Calculations 

presented assume that the nozzle was unconstrained from expansion, and that the annulus and shell imposed no 
stresses on the nozzle, having assumed that the graphite gasket is able to absorb them. 
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Figure 6. CFD calculated surface heat transfer coefficient 
plotted as a function of axial distance along the burner. 


Figure 7 Finite elements solution for temperature distribution in burner assembly 


One representative calculation is presented here (Fig. 7) in which the shell, nozzle and annulus are in the steady- 
state condition. The flow surface heat transfer is modeled using Equation (1). The surfaces of the water-cooling 
passages are set to 373 K, while the annulus radiates to the surroundings, assumed black at 300 K. Heat transfer 
between the nozzle and the annulus is by both conduction and radiation. Emissivity of the nozzle is assumed 0.9 
and of the annulus 0.7. Components conduct where they contact, except that a thermal resistance 0.0001 m 2 K/W is 
placed between the nozzle and the annulus and between the nozzle and the shell to model a graphite gasket of 0.5 
mm thickness. No flow is assumed in the coflow passage. This latter point is important since, in all our experiments 
described below, there is a coflow, which cools the nozzle and annulus. 

The maximum temperature in the steel is 1180 K, below the maximum service temperature in air (which is 1310 
K for intermittent service), while the maximum temperature in the nozzle is reached at the exit of the nozzle and is 
1680K, below the maximum for the silicon carbide (which is -1900 K in an oxidizing atmosphere). In this 
condition, the maximum positive value of the first principal stress in the nozzle occurs in the outer part of the flange 
and is -30 kpsi (206 MPa). The tensile strength of silicon carbide is represented in the published data as a 
cumulative probability of failure. 14 Although there is scatter in the data and a trend of increasing strength with 
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temperature, the probability of failure typically rises above 0 at ~22 kpsi (150 MPa)and approaches 1 at ~44 kpsi 
(300 MPa). Thus there is a chance of nozzle failure under the calculated conditions. High thermal stresses are 
associated with high temperature gradients, particularly in the radial direction. Note that with flow in the coflow 
nozzle the nozzle temperature is reduced to -1200K -1400 K, reducing temperature gradients and thermal stresses to 
acceptable levels. Increasing the contact thermal resistance of the gasket increases the temperature of the nozzle in 
the flange area, but also reduces temperature gradients and thermal stresses. 

Time-dependent calculations were performed in which the temperature of the assembly was initially 300 K. 
These showed that, while the component temperatures increase monotonically towards the steady state, the 
maximum stress peaked at as much as two times the steady state value at around 10 s, before falling to steady state. 
This is because thermal stresses depend upon temperature differences within the nozzle, especially in the radial 
direction, and upon the thermal expansion coefficient. The nozzle flow surface temperature rises more rapidly than 
temperatures some distance from the surface, causing an overshoot in temperature differences. The problem is 
worsened by the increase in thermal expansion coefficient with temperature. 

To avoid failure in the nozzle, the burner, in practice is always operated with coflow. No problems with the 
nozzle have been encountered over a wide range of conditions. If operation without coflow were to be attempted at 
the high-heat-flux conditions then the graphite gaskets would have to be replaced with insulating gaskets from a 
material such as mica. For operation without coflow it would be desirable to water cool the annulus. 
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Figure 8. Comparison of attached flames as 
observed in the infrared: (a) base attached 
flame (M= 1, 7M700K, vitiated H 2 with air 
coflow) and (b) external attached flame (M= 1, 
M ent hai py = 7, H 2 vitiated air with H 2 coflow). 


Figure 9. Matrix of burner operating conditions 
with several nozzles (M= 1, 1.6 and 2). H 2 -vitiated 
air center jet with H 2 coflow. Infrared (false color) 
images. Red box refers to case with CARS data. 


IV. Flow Visualization of Test Matrix 

The flowfield at the exit of the burner was visualized by employing either an infrared (IR) camera or a visible 
light color consumer-grade digital camera (and sometimes both). The IR camera was a FLIR Systems Thermacam 
SC3000 imaging in the 8 - 9 micron (long wave) region; the effective shuttering time of this camera is significantly 
less than 0.3 ms, allowing some of the unsteady structure to be resolved. The visible light camera exposure time 
was - lsec, long in relation to the unsteady structure. Intermittent operational problems with the IR camera required 
that the visible camera be used as a second choice on some runs. 

Figures 8-12 contain tables and flow images showing the various types of flames observed in the test matrix. 

The tables contain information pertaining to the following states of the flame: no flame (NF), detached flame (DF), 
flame attached at the base (BAF) or at the external coflow boundary (EAF), and whether the flame holding is 
marginal, i.e., at the point of extinction (M). Figure 8 contains IR images that illustrate the two types of attached 
flame. The BAF in Fig. 8(a) is a flame that extends to the annular base region between the central and coflow 
nozzle exits, where the flow is recirculating. The existence of a flame here indicates that the ratio of fuel to air that 
enters this region through mixing is within the range of flammability. This type of flame seems only to be observed 
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with vitiated hydrogen in the center jet. The EAF appears to be a diffusion flame that reaches to the burner surface 
at the boundary between the coflow and the ambient air. This flame is possible only in cases where the coflow is a 
fuel, occurring because of the low subsonic speed of the coflow. Notice some asymmetry in the EAF flame of Fig 
8(a): this is believed to be due to the burner being assembled with the silicon carbide nozzle slightly off center, 
causing the width of the coflow nozzle exit to vary slightly around the circumference. 



M= 1 

M= 1.6 

M= 2 

M enthalpy 

ER=1 

ER=1 

ER=1 

5.5 

6 

6.5 

NF 

DF 


NF 

NF 

NF 

7 

DF 


uf 

NF 



M = 1 

M=1.6 

M= 2 

M enthalpy 

ER=1 

ER=1 

ER=1 

5.5 
6 

6.5 
7 

0 

EAF^ 


0 

EAF 

/EAF 

Ieaf 



33% H 


Figure 11. Matrix of burner operating 
conditions with C 2 H 4 vitiated air center jet 
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Figure 9 shows IR images for cases with H 2 - vitiated air in the center jet and H 2 in the coflow. The sensible 
enthalpy of the center jet is varied from that of 
Mach 5 flight to that of Mach 7; the exit Mach 
number of the center jet is 1.0, 1.6, and 2.0; and 
the equivalence ratio is 0.5 and 1. For M= 1 
and 1.6 an exterior attached flame is observed 
in all cases, whereas for M = 2 there is either no 
flame or the flame is detached (stands off from 
the nozzle exit). Where the flame is detached it 
spreads rapidly, perhaps being coupled with 
shock waves. The trend with increasing center 
jet enthalpy for a detached flame is for the 
flame to move towards the nozzle exit. 

Figure 10 shows both IR and visible images 
for cases with H 2 -vitiated air in the center jet 
and C 2 H 4 in the coflow. The IR images 
indicate a detached flame for cases at center jet 
Mach numbers of 1 and 1.6, and at the highest 
enthalpy. The visible images show that the 
detached flame is associated with a deep blue 
emission, perhaps CH or C 2 emission. In other 
cases there is no combustion in the jet, although 
some emission is observed, indicating perhaps 
some dissociation of the C 2 H 4 . Comparison of 
the table entries with Figure 9 shows that C 2 H 4 
has poorer flameholding than H 2 , as is expected Figure 12. Matrix of burner operating conditions with 

based on measurements of extinction. 20 vitiated H 2 center jet and sonic air coflow. Visible light 

images. 


900 K 
1000 K 
1100 K 
1200 K 
1300 K 
1400 K 
1500 K 
1600 K 
1700 K 
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Figure 1 1 shows both IR and visible images for cases with C 2 H 4 -vitiated air in the center jet and H 2 in the 
coflow. The IR images indicate an exterior attached flame for all cases that lighted. The limiting factor is whether 
or not combustion takes place in the burner. The 0’s in the matrix indicated that it would not take place. Where it 
was possible to light the burner at all it was lit at the Mach 7 flight enthalpy condition; for the M= 2 nozzle case, 
once lit, it was found to be possible to reduce the enthalpy and the burner remain lit. Emissions levels in the center 
jet are higher than for the H 2 -vitiated cases. In the IR this may reflect emission from the C0 2 . In the visible this 
may reflect the presence of soot or soot precursers. Comparison with Figure 9 shows that H 2 has better flame 
holding with the C 2 H 4 -vitiated air than the cases with H 2 vitiation. This may be due to higher concentration of 
radical species or soot precursors in the center jet. If true, the results with this small burner may not be 
representative of the effect of vitiation in large scale hypersonic propulsion test facilities, where, presumably, the 
facility bums the fuel to a nearly equilibrium state. 

Figure 12 shows visible light images for cases with vitiated H 2 in the center jet, either 33% or 50% by volume of 
the jet flow being H 2 . The flame is made visible by a pale blue emission from the OH radical. In some cases, no 
flame is observed, in others detached flames or flames attached at the base. The trend with increasing T is from no 
flame, to a detached flame, and then to a base-attached flame. There seems to be no consistent effect of M on flame 
holding, although at higher M the flame is larger downstream due to the greater flow rates of both vitiated air and 
coflow H 2 . One feature of these flames was a high level of noise. This noise increased as the enthalpy was reduced 
and the flame became more unstable. Also, the quality of the noise changed, from white noise for an attached flame, 
to a low pitched “rumbling” sound as the flame became detached, and finally to a “crackling” noise before the flame 
went out. This severe noise is neither conducive to operator nor to optical system presence in the laboratory. 

V. Optical Diagnostics Measurements 

A. The Dual-Pump CARS-IRS Method 

The dual-pump CARS system 21 uses a seeded Nd:YAG laser at 532 nm (green) and a narrow-band dye laser at 
553 nm (yellow) as pump beams and a spectrally broad dye laser centered at 607 nm (red) as the Stokes beam. 

Some of the green laser is used for pumping the yellow pump and red Stokes dye lasers. The frequency difference 
between the green and red beams corresponds to the vibrational Raman shift of N 2 . The wavelength of the yellow 
pump beam is such the frequency difference between the yellow and red beams equals the vibrational Raman shift 
of 0 2 . Therefore the CARS system probes N 2 , 0 2 , and also H 2 because pure rotational H 2 transitions are present in 
both spectral regions. The laser beams are combined at the focusing point of a spherical lens in a folded BOXCARS 
geometry, generating a measurement volume of about 1.5 mm x 0.2 mm x 0.2 mm. The coherent CARS signal beam 
generated at 49 1 nm (blue) is collected and collimated by another spherical lens, focused on the input slit of a one- 
meter spectrometer, and the spectra are recorded on a cooled CCD camera. Post-processing software fit spectra with 
a theoretical model. The shape of these spectra provides information on the temperature while the relative intensities 
of these spectra provide a measure of the relative mole fractions. 

The interferometric Rayleigh scattering system 9,10,21 uses as the light source the same green laser used by the 
CARS system. To increase the Rayleigh signal and to avoid laser- induced gas breakdown, an energy storage delay 
line of about 1 8 nsec (pulse stretcher) is used in the green beam optical path. The receiving optics collect Rayleigh 
scattered light from two directions, the scattered light from the yellow and the red CARS beams is removed, and the 
remaining green light is passed through a solid etalon. To provide a reference, green light directly from the laser 
(not Rayleigh scattered) is also injected into the etalon. An electron multiplication gain CCD camera records the 
interference fringe pattern, which is formed at the output of the etalon. Interferogram processing software extract 
Rayleigh spectra and fitting algorithms solve for Doppler shift and implicitly for velocity. 

B. Preliminary Results 

Figure 13 shows some preliminary dual-pump CARS measurements acquired along a diameter at 2 mm above 
the nozzle exit, with the M=1.6 nozzle, H 2 -vitiated air at M entha i py =6 and H 2 coflow. A visible light image of this 
flame is inset in Fig. 13(a) - see also the matrix in Fig. 9. The nozzle exit conditions calculated with ID flow 
assumptions, and assuming 20% heat loss (see Section IIIB) are: V= 1064 m/s, T= 1094 K, and mole fractions of 
N 2 and 0 2 of 0.589 and 0.193 respectively. For this data set, the CARS measurement volume was scanned at a 
constant rate across the jet diameter in a direction parallel to the laser beams while acquiring 1000 points at 20 Hz. 
Single shot (instantaneous measurements) are plotted. In the nozzle base and coflow jet regions, in the range 5 mm 
< | x | < 9 mm, there was insufficient N 2 to obtain meaningful fits of the CARS spectral to obtain temperature and 
composition. No data is shown in these regions. This absence of N 2 (and also 0 2 ) in the base region explains the 
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lack of any flame here. (As previously noted, this case had an external attached flame.) Elsewhere, 94% of the data 
were successfully fit. The data are seen to scatter around the trend lines. This scatter is due in part to instrument 
random error, and in part to turbulent fluctuations. The temperature varies symmetrically, from about 1400 K at the 
axis down to less than 800 K at the edge (x = ± 5 mm). This variation is attributed to heat transfer to the burner and 
nozzle walls. The ID static temperature based on an assumed heat loss of 20% is 1094 K; if the assumed heat loss is 
reduced to zero, the ID static temperature is 1300 K, still 100 K below the centerline observation. On the other 
hand, the mole fractions of 0 2 and N 2 are fairly flat in the center (x = ± 4 mm) region of the jet, and compare well 
with the ID values. At the outer edge of the coflow jet (x = 9 mm) there are some high temperatures, which may be 
the external attached flame. However, the data here should be treated with suspicion since the spatial resolution is 
insufficient to resolve this part of the flow. 



Figure 13 Dual-pump CARS measurements along a diameter at 2 mm above the nozzle exit, with M= 1.6 H 2 - 
vitiated air at M entha i py = 6 and H 2 coflow, (a) Temperature (inset is visible light image of flame), (b) species 
mole fractions. 


Interferometric Rayleigh scattering data were not acquired at 2 mm above the flow due to interference of 
scattered light from the nozzle exit. Data were acquired (simultaneously with CARS) at other locations. This data is 
still being analyzed, although some is presented in the companion paper by Tedder et al. 21 . The intent herein is to 
present data that can be compared with calculations for the purpose of verifying our techniques. 


VI. Summary and Conclusions 

This paper describes progress towards generating an experimental data set that can be used by CFD modelers to 
develop models for turbulent supersonic combusting flows. A laboratory supersonic scale burner is described which 
duplicates, in essential features, the flow of a larger facility in which the final experimental data is to be acquired. 
The flowfield consist of an axisymmetric free jet of hot reaction product and an unheated annular coflow jet. By 
selection of the composition of burner reactants and coflow, a supersonic free jet with combustion can be created. 
This burner has been designed to verify the choice of experimental test conditions of the final experiments, and to 
verify the optical diagnostic instruments that will be employed. The calculated range of available flow conditions of 
the burner is described, along with some detailed design calculations. The burner is operated over a wide matrix of 
available conditions, and the supersonic combusting jets are visualized in the visible and infrared wavelengths to 
assess flame holding. The final experiment flow field will be mapped using the dual pump CARS-IRS technique, 
which can provide simultaneous temperature, composition, and velocity. From series of such measurements, many 
of the turbulence statistics of interest to modelers can be found. The fluid dynamical length scales of the laboratory 
facility are not small in relation to the CARS system measurement volume; hence the scale of the final experiment 
will be more than six times larger. However, the small burner may provide useful data using other techniques with 
better spatial resolution. In this paper we have discussed the CARS-IRS technique and presented a sample of some 
preliminary experimental data acquired in the laboratory burner to verify it. 

The ability to create a supersonic combusting jet suitable for study with optical techniques has been 
demonstrated. These cases have fully compressible mixing layers. Cases with flame holding and with flame 
standoff, as well as with both H 2 and C 2 H 6 fuels are available. The dual-pump CARS-IRS technique has been 
demonstrated, and some shortcomings identified. 
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